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Abstract A detailed account is given on approximation schemes to the Ein- 
stein theory of general relativity where the iteration starts from the Newton 
theory of gravity. Two different coordinate conditions are used to represent 
the Einstein field equations, the generalized isotropic ones of the canoni- 
cal formalism of Arnowitt, Deser, and Misner and the harmonic ones of the 
Lorentz-covariant Fock-de Donder approach. Conserved quantities of isolated 
systems are identified and the Poincare algebra is introduced. Post-Newtonian 
expansions are performed in the near and far (radiation) zones. The natu- 
ral fitting of multipole expansions to post-Newtonian schemes is emphasized. 
The treated matter models are ideal fluids, pure point masses, and point 
masses with spin and mass-quadrupole moments modelling rotating black 
holes. Various Hamiltonians of spinning binaries are presented in explicit 
forms to higher post-Newtonian orders. The delicate use of black holes in 
post-Newtonian expansion calculations and of the Dirac delta function in 
general relativity find discussions. 



1 Introduction 

In the weak-field slow-motion limit of any theory of gravity the Newtonian 
theory of gravity comes into play because it describes the motion and struc- 
ture of gravitating objects very well in that regime. According to current 
knowledge from experiments and observations, the most reliable theory of 
gravity is the Einstein theory of general relativity, [T] . 
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Most of the objects in the Universe seem to have velocities v which are 
small compared with the speed of light (in the following denoted by c), i.e. 
v/c < 1/3. In those cases it is to be expected that the Newtonian theory is 
a good starting point for an iteration scheme toward full general relativity 
thinking in terms of expansion of general relativity in powers of dimensionless 
v/c. The conservative nature of most of the phenomena, i.e. no measurable 
gravitational radiation damping, the power series has to be a one with or- 
dering parameter (v/c) 2 (recall v 2 /c 2 < 1/10 in the regime of our interest) 
because of motion-inversion symmetry in those cases. We shall see later that 
at the order (v/c) 5 the gravitational radiation damping enters for the first 
time. Calling the order (v/c) 2n the nth post-Newtonian order (n PN), the 
gravitational dissipation (because of gravitational radiation damping from 
radiation emission) enters at the 2.5 PN order so the exponent n will also in- 
clude half-integer numbers. In case of bound systems the virial theorem holds 
which tells us that v 2 ~ GM/r, where G is the Newtonian gravitational con- 
stant, r a typical distance of two bodies or a typical radius of one body, and 
M the total mass of the system. So an expansion of general relativity in pow- 
ers of (v/c) 2 is at the same time an expansion in powers of GM/c 2 r if bound 
systems are considered. From the dynamics of bound systems the dynamics 
of low-velocity scattering is straightforwardly obtained. 

Depending on the mathematical representation of the Einstein theory 
(choice of coordinates, choice of variables, etc.) and on the physical aspects 
under investigation (physics in the near zone, physics in the far zone, etc.) 
there exist many of those expansions which all are called PN expansions. The 
crucial ordering parameter in all PN expansions is always 1/c 2 . The Einstein 
field equations as well as the equations of motion which follow from them 
by integrability conditions can be formally expanded in powers of 1/c 2 . The 
solutions of these equations, however, can not be expanded in this way, in 
general. Only under conditions where no gravitational radiation is present 
an all-over-in-space expansion in powers of 1/c 2 is feasible. This is the case 
for stationary systems only. Yet for those parts of a radiating system where 
radiation plays no role, i.e. in the conservative parts, an expansion in powers 
of 1/c 2 can be carried out too. In general, PN expansions in powers of 1/c 
are feasible in the near zone (r < < A with A a typical wavelength of the grav- 
itational radiation and r the radius of a sphere enclosing the matter source) 
and in the far or radiation zone (r » A) of a radiating system but they are 
not analytic in 1/c because c is showing up in log-terms. In the near zone, 
up to 3.5 PN order, which means order 1/c 7 = (1/c 2 ) 3 5 , a PN expansion in 
powers of 1/c is valid, where the even powers are connected with the conser- 
vative dynamics and the odd powers with the dissipative radiation reaction 
dynamics. From the 4 PN order on (corresponding to 1.5 PN order in the far 
zone), log-c terms enter via tail effects which emerge from back scattering, 
[2]. Post-Minkowskian (PM) series are global series in space where only the 
weak-field limit is assumed without restrictions to the velocities of the bodies. 
Here, the ordering parameter is G or, dimensionless GM/c 2 r, which is equiv- 
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alent to a non-linearity expansion of the Einstein field equations, [2\. If the 
virial theorem holds, namely (GM/c 2 r) ~ (v/c) 2 , PM series can be further 
expanded into PN series. Having much less information on their respective 
orders, the PN series are much easier to be worked out analytically, compare 
e.g., the 1 PN point-mass Hamiltonian, see e.g. [3], with the 1 PM point-mass 
Hamiltonian, only recently achieved in [4]. 

The most extensively performed PN calculations in the literature can be 
divided into two classes, the class applying harmonic coordinate conditions, 
see e.g. [S], and the other class using Arnowitt-Deser-Misner (ADM) gen- 
eralized isotropic coordinate conditions [6]. Whenever detailed comparisons 
between the two classes have been made, agreement of the results could be 
achieved, see e.g. [7], [S]- This is quite an important aspect in view of the 
use of Dirac delta functions for modelling black holes in explicit calculations. 
These functions are a-priori problematic in a non-linear theory like general 
relativity so sophisticated regularization methods have to be employed. The 
only regularization method which turned out to be successful to the highest 
orders of the explicit calculations is the dimensional one, [9], [10] . In sect. 13.41 
a detailed account of black holes together with their representing Dirac delta 
functions will be given because they blatantly violate the weak-field condi- 
tion. In this case, only external velocities of the black holes can be small 
and their relative gravitational interaction weak, and only in this sense black 
holes and Dirac delta functions can fit into PN schemes. 

The binary point-mass dynamics at 3 PN order and the gravitational waves 
emitted hereof (3 PN wave generation) could have been managed successfully 
through only dimensional regularization, [S], [TU], [IT]. Extended body cal- 
culations have been performed too but they are much more complicated and 
still not fully under control at 3 PN order [12] . Only a surface- integral-based 
method of the Einstein-Infeld-Hoffmann type did succeed too, [13], [14] . 

The many-body systems treated most successfully with PN approximation 
methods arc point-mass systems composed from two or more objects and 
spinning binary systems with even rotationally deformed components; for 
the most recent results with spin, see [TS], [IB]. [T7j. [I8j. [i"9j. [20j, [2T], 

Recently, an effective field theory approach has been advocated for PN 
calculations [21]. Applications to point-mass systems and to systems with 
spinning components have already been performed through higher PN orders, 
[25] , [2~4] , [25] . An obvious difference in the spin dynamics between [2T] and 
[24] found clarification in [26] showing the correctness of the result in [21] . 

In this article, Latin indices from the mid alphabet are mostly running 
from 1 to 3, v = (i/),v 2 = v i v % , dt and di denote the partial derivatives 
with respect to the time and space coordinates t and x l , respectively, and 
the functional derivative (Frechet derivative) of a functional, say F[f], with 
respect to a function / takes the form 5F[f]/Sf. Greek indices are mostly 
running from to 3, whereby x° — ct. The signature of the four-dimensional 
metric is +2. Particles arc numbered with indices from the beginning of 
the Latin alphabet. 
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2 Systems in Newtonian gravity in canonical form 



In the Newtonian theory, the equations governing the motion of gravitating 
ideal fluids are (i) the equation for the conservation of mass, 

d t Q* + div(g*v) = 0, (1) 

where is the mass density and v the velocity fied of the fluid, (ii) the 
equations of motion, 

£*<9tv + ^ grad v 2 - v x curl v = -grad p + grad U , (2) 

where p is the pressure in the fluid and U the gravitational potential, and (iii) 
the equation of state, using internal energy density e and specific enthalpy h, 

e = e(f>*, s) with de = hdg* + g^Tds, or dp = g*dh — g*Tds , (3) 

with (iv) the conservation law for the specific entropy s along the flow lines, 

d t s + v • grad s = . (4) 

The gravitational potential reads, 

tf(x,t)=G /'rf 3 x'^4r. (5) 
J |x x | 

It results from (v) the Newtonian field equation, 

AU = -AitGg* . (6) 

A is the Laplacian and |...| means the standard euclidean distance. 

Written in form of Hamilton equations of motion, i.e. d t A(x,t) = {A(x,t),H}, 
the equations above take the forms, (i) the mass conservation equation 

di j-g* , (7) 



dt \ 8-Ki 

notice v l = (ii) the equations of motion 



d ^ a ( 5H \ « ( 5H \ « ( SH \ , SH 



and (iv) the entropy conservation law 
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where the Hamiltonian is given by H — H(g*, 7T,, s) with 717 the linear momen- 
tum density of the fluid, see [57]. Here, use has been made of the kinematical 
Lie-Poisson bracket relations between the fundamental variables 

{m (x, i) , g« (x' , t ) } = ^- [g* (x' , t)S (x - x' )] , (10) 

{^(x, t), S (x', t)} = ^r-si.* - x') , (11) 

{^(x, t), Kj (x', t)} = tt^x', *)^j«(x - x') - vr^x, t)^-8{^ - x') , (12) 

and zero otherwise, where #(x— x') denotes the standard Dirac delta function 
in three-dimensional space. It fulfills J d 3 x <5(x) = 1. 

In the Newtonian theory, the Hamiltonian of the fluid is given by, 

If f W _ G f d 3 xrf 3 x ,Mx,t)Mx',t) /■ ^ 

27 g* 2 J |x-x'| y v y 

For point masses, the total momentum and mass densities read, consistent 
with the Eqs. OH) and <H2J), 

7ij = y^p g j^(x-x ct ) i g* =y^m g ^(x-x a ), (14) 

a a 

where the position and momentum variables fulfill the standard Poisson 
bracket relations, 

{x l a ,p a j} = Sij, zero otherwise , (15) 
and the Hamiltonian takes the form, 

2 V rna 2^|x a -x b |' 

where the self-energy term has been dropped (for regularization techniques, 
see section! 



3 Canonical general relativity and PN expansions 

In curved spacetime the stress-energy tensor of an ideal fluid takes the form 

T» v = g{c 2 + h)u»u v + pg^, g^u v = -1 , (17) 

where g denotes the proper rest-mass density, h the specific enthalpy, and 
it^ the four-velocity field of the fluid. Using energy density e = g(c 2 + h) — p 
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(also the specific internal energy 77 = ej g — c 2 could be used), the equation 
of state reads 

e = e(g, s) with de — (c 2 + h)dg + gTds, or dp = gdh — gTds . (18) 
The variables of the canonical formalism are chosen to be 



gu°g, s, 



1 



7r t = -J—gTl (19) 



c 



They fulfill the same (universal) kinematical Lie-Poisson bracket relations as 
in the Newtonian theory, see [27], or also [28] . 

{7Ti(x, <), e*(x', *)} = ^-b*(x', t)<J(x - x')] , (20) 

{ 7ri (x,t), S (x',t)} = ^^«5(x-x'), (21) 

'^(x-x')-^(x,t)^. 
The evolution equations take the form 



{TTifotJ.TTjCx'.t)} = ^(x',t)— <5(x - x') - ^(x,i)— «J(x - x') . (22) 



^ = -9i(^) = <VY~i'«">=(>. (2:il 



^=-^9 iS <=► u^ fl = 0, (24) 

OT 07T ; 



corresponding to <9 M (^Tf) - - y^gT^ dig^ = , 

v = - — , where v = c—z . (26) 

07Tj U u 

The linear and angular momenta of the fluid read, respectively, 

Pi = J d 3 x 7Ti, J* = J d 3 x e ijk x j TT k . (27) 

For a system made of point masses simplifications take place, 

h = p = s = , (dusty matter), (28) 

and further, 
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£>* =^TO a (5(x-X a ), TTi=^2pai6(-K-X a ), v\ = -r± , (29) 

a a 

where p a i and a;* respectively are the linear momentum and the position 
vector of the ath particle. The kinematical Poisson bracket relations are given 
by 

{x l a ,p a j} = Sij, zero otherwise . (30) 
Hereof the standard Hamilton equations result, 

dPaA = _dH_ dxi = dH_ 
dt dx\ ' dt dpai ' 

Remarkably, the difference to the Newtonian theory comes solely from the 
Hamiltonian which is thus a dynamical difference and not a kinematical one. 
This statement refers to the matter only and not to the gravitational field. 
The latter is quite different in general relativity. 



3. 1 Canonical variables of the gravitational field 

Within the ADM canonical formalism of general relativity, in generalized 
isotropic coordinates, the independent gravitational field variables hj? and 
7r^ T enter in the form 

9H = ( 1 + I A + hj? , (32) 



7r y ' = n ij + 7T^ T , (33) 

where gij — gji = jij is the metric of the curved three-dimensional hy- 
persurfaces t = const, it^ c 3 /16itG is the canonical conjugate to 7^, i.e. 
7r y = — ^l^iK^ — ^K^), where Kij — Kji is the extrinsic curvature of 
the t — const slices, and 7 = det(Tjj), 7**7y = ^*j> an d f° r ^ holds 

f« = d i7 r j + dy - ^ 3 d k ir k . (34) 

Obviously ir u = 0, or n\ — n^hf?. The canonical conjugate to hj? reads 
7r^ T c 3 /WwG. The index TT means tranverse-traceless, i.e. hj^ = n^ T = 0, 
djhj^ = djir l 4 T = 0. 

Using those variables, the Einstein field equation y/^gG 00 = ^r-y/^gT 00 
can be put into the form, employing point masses for the source, 
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1 / 2 R = 



7 



7T*7r? 7r?7T? 

1/2 \ 3 I 2 3 



16ttG 



X! ( TO a c2 + l l3 PaiPaj) 1/2 S a ,(35) 



and the field equations ^gG\ = ^f-y/^T° read 

16ttG 



2d ] TT 3 +7T kl d ll kl = 



^Pai&a , 



(36) 



where S a = 5(x — x a ). The Eqs. (|3"B"|) and (f3"6")) are the famous four constraint 
equations of general relativity. 

In the gauge Eqs. (|32|) - (|34|) the ADM Hamiltonian can be written, [6], 



H 



x ai Pat i "-ji ) " ' 



ij ' "TT 



16ttG 



(37) 



resulting from the solution of the four (elliptic-type) constraint equations. 
The additional Hamilton equations of motion for the gravitational field are 
given by 



dt 



16ttG SH 

C 3 JftTT' 



Oh 



TT 



16ttG <5# 
c 3 j7r^ T 



(38) 



The transition to a Routh functional simplifies a lot the construction of the 
dynamics of the matter and of the gravitational field. The Routh functional 
is chosen in the form, [29] . 



R[xi,p ai ,hJ?,d t hJ?] =H 



16ttG 



d 3 x 7T^ T 9 t /l^ T . 



(39) 



The evolution equations for the matter and the gravitational field now read 



5jR(t')dt' 
5hJ?(x k ,t) 



Pai 



dR 

dxl 



dR 

Opai 



(40) 



The conservative dynamics results from the on-field-shcll Routh functional 



with solved field equations, in the form 

. 6 J RebenWdt' . u 8 J R shen (t')dt' 

Pai{) ~ Sxi(t) ' XAt) - 5p ai (t) 



where 



S J Rshel\(t')dt' 9i? sh oll d <9-R sh ell 



5z{t) 



dz{t) dt dz{t) 



+ ..., z=(x l al p ai ) 



(41) 



(42) 



(43) 
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Using the matter equations of motion in the Routhian i? s heii the Routhian 
can be brought into the form R(x l a ,p a i). Herein, however, the meaning of the 
variables x l a and p ai has changed, see [30], [3T], [52] . 



3.2 Brill- Lindquist initial-value solution for binary 
black holes 

The Brill-Lindquist solution for multiple black holes is a pure vacuum solu- 
tion of the constraint equations at initial time t under the conditions of time 
symmetry, i.e. p a i = = ir^ , and of conformal flatness, i.e. hj? = 0, [33] . 
A related vacuum solution is the Misner-Lindquist solution where an addi- 
tional isometry condition is imposed, [34] . [35] , Under those conditions, the 
only remaining constraint equation reads, not using vacuum but (point-mass) 
sources, 

1 ^^=^r mA , ^7 = = (44) 



S 



In the case of two black holes, its solution is given by, see [36] . 

AG ( OL\ Qf2 

' — + — 

T\ r 2 



(45) 



with (a, 6=1,2 and b ^ a) 

"»= o + h/ 1 + -a- ,~ + o-a- ,~ - 1 i ( 46 ) 




resulting into the Brill-Lindquist solution for binary black holes. Obviously, 
each Brill-Lindquist black hole is represented by a Dirac delta function (ficti- 
tious image mass-point; see section [3~4]) . In the Misner-Lindquist case, infinite 
many fictitious image mass-points are needed for each black hole, [34] . [35] . 
[36]. 

The energy of the Brill-Lindquist solution simply reads 

#bl = (ai + a 2 )c 2 = (mi + m 2 )c 2 - G , (47) 

The methods which have been used for the obtention of the Brill-Lindquist 
solution from sources (notice, in the original work of Brill and Lindquist this 
solution has been obtained without any regularization as a purely vacuum 
solution) are analytical Hadamard regularization and mass renormalization, 
[36] , as well as dimensional regularization based on the d- dimensional metric 
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/ 1 d — 2 \ 3=5 
with solution [T denotes the Euler gamma function) 

(for more details see section l3"Hj) . The PN expansion of the Brill-Lindquist 
initial energy expression is straightforward to all orders of 1/c 2 . Once it has 
fixed the static ambiguity parameter u static (see [32) ) in non-dimensional- 
regularization calculations to the correct value of zero, [36]. At that time, 
however, it was not quite clear that the Brill-Lindquist solution delivers the 
correct boundary conditions for the point-mass model. 

The truncation of the constraint equations in the form hj^ = as well as 
dropping an additional term in the Hamiltonian constraint connected with 
the energy density of the field momentum results in a remarkable, fully ex- 
plicitly solvable conservative so-called skeleton dynamics which allows a PN 
expansion of the Hamiltonian, and of all the metric coefficients too, to all 
orders, see next section. 



3.3 Skeleton Hamiltonian 

In Ref. [37] the skeleton dynamics has been developed. The skeleton approach 
to general relativity requires the conformal flat condition for the spatial three- 
metric for all times (not only initially as for the Brill-Lindquist solution) 

Hj = (1 + I0) 4 % • (50) 
Hereof, in our coordinate system, maximal slicing follows, 

7^=2^77^ = 0. (51) 

Under the conformal flat condition for the spatial three-metric, the momen- 
tum constraint equations become 

a 

The solution of these equations is constructed under the condition that 7rf 
(and not 7r y , see Eqs. ([33| and (f34|) l is purely longitudinal, i.e. 



(53) 
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This condition is part of the definition of the skeleton model. At spacelike 
infinity, the surface-area integrals of irj or tt ij are proportional to the total 
linear momentum of the binary system. 

Furthermore, in the Hamilton constraint equation, which in our case reads 

Arh= 16^?^ mj a f v\ \ 1/2 

* (1 + W * (1 + i 0) + l^mlc- ) ' [b } 

we perform a truncation of the numerator of the first term in the following 
way 

= -Wjdrf + fl((2^7rj) -» -IV J)-' = ^l^PajVjSa , (55) 



i.e. we drop from 7r?7r*- the term di(2VjTTj). This is the second crucial trunca- 
tion condition additional to the conformal flat one. Without this truncation 
neither an explicit solution can be achieved nor a PN expansion is feasible. 
From we know that at the 3 PN level the /i^ r -field is needed to make 
the sum of the corresponding terms from irf 7rJ analytic in 1/c. 
With the aid of the ansatz 

a a 

and by making use of dimensional regularization, the energy and momentum 
constraint equations result in an algebraic equation of the form, [37] , 



1 + 

r a b 



Pl/(m 2 a c 2 ) 



(1 + Aa b /r ab f 



PaiVai/c (57) 



{l + Aa b /r ab ) T 



where A = G/(2 c 2 ) and b ^ a. 

With these inputs the skeleton Hamiltonian for binary black holes becomes 
(at least initially, for p a — 0, the solution is consistent with general relativity) 



16ttG 



J d 3 xA(f> = c 2 ^a a . (58) 



The Hamilton equations of motion read 

Xa = ^' Pa = -thT a - (59) 

In the center-of-mass frame of the binary system, we define 

p = pi=-p2, TEEEXi-x 2 , r 2 = (xi - x 2 ) ■ (xi - x 2 ) . (60) 
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Further, we will employ the following convenient dimensionless quantities 
tc 3 „ rc 2 , p - i? s k 



# s k = -5 , (61) 
... /*r 

J c „ p r 



G m G to /ic /1 c z 



i = 7f— , P, = ^, P 2 =P 2 +) 2 /f 2 , (62) 

where J = r x p is the orbital angular momentum in the center-of-mass frame 
and p r = p • r/r the radial momentum. The total rest-mass is denoted by 
m = mi + m,2 and the reduced mass by /1 = mima/m. The binary skeleton 
Hamiltonian Hsk can be put into the following form, |38j , 

ff sk = 2f( V>i + ^2 - 2 ) with (63) 



4 ^ 2 (p 2 . + j 2 /f 2 ) \ 1/2 (8p 2 + 7j 2 /f 2 ^) i^ 2 



^i = 1 + t4-(i + Srn "J «~2/7 ' >( 64 ) 

4 r ip2 \ X- V2 / ° ^ ^2 



x+ / 4^ 2 (p 2 +?/f 2 )x 1/2 (8j5 2 + 7j 2 /f 2 ) !/ 2 
4fVi V 1+ X+^i J 8f 2 v>I 



where \- = (l — \A — 4 1/) and x+ = (l + V 1 — 4 1/) with f = /j/to. 

The conservative skeleton Hamiltonian has the following nice properties. 
It is exact in the test-body limit where it describes the motion of a test 
particle in the Schwarzschild spacetime. It is identical to the 1 PN accurate 
Hamiltonian for the binary dynamics in general relativity. Further, as ex- 
plained earlier, when point particles are at rest, the Brill-Lindquist initial 
value solution is reproduced. It is remarkable that the skeleton Hamiltonian 
allows a PN expansion in powers of 1/c 2 to arbitrary orders. The skeleton 
Hamiltonian thus describes the evolution of a kind of black holes under both 
conformal flat conditions for the three-metric and analyticity conditions in 
1/c 2 for the Hamiltonian. Of course, gravitational radiation emission is not 
included. It can, however, be added to some reasonable extent, see [38] . 

Restricting to circular orbits and defining x = (Gtow/c 3 ) 2 / 3 , where u is 
the orbital angular frequency, the skeleton Hamiltonian reads explicitly to 3 
PN order, 

x /3 v \ 2 /27 29 17 ,\ , 

H sk = hi — I ar + 1 v v 1 \x i 

2 \8 24 J \16 16 48 J 

675 8585 7985 , 1115 ,\ , m . 5 . 

H v v H v 3 ]x i + 0(x 5 ). (66) 

128 384 192 10368 J 

In Ref. [37] the coefficients are given to the order x 11 inclusively. 
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In the Isenberg- Wilson-Mathews approach to general relativity only the 
conformal flat condition is employed. Thus the energy stops being analytic 
in 1/c at 3 PN. Through 2 PN order, the Isenberg- Wilson-Mathews energy 
of a binary is given by 

x /3 v \ 2 (H 39 17 ,\ , 
FlWM = 2 + U + 2i> + U - iT ~ 48* ■ (67) 

We already quote here the 3 PN result of general relativity. It reads, see Eq. 

x ("S v\ 2 (21 19 1 ,\ 3 

H-iPN = M-H ar + ^ H ^ *c 

J 2 V8 24/ V 16 16 48 

205 2 34445 \ 155 2 35 .. , , 




192 



it 2 \v+ v 2 + v A he 4 . (68) 

1152 J 192 10368 / v ' 



The difference between Hiwu an d H s ^ through 2 PN order shows the effect 
of truncation in the field-momentum part of i? s k and the difference between 
ifiWM an d Hzpn reveals the effect of conformal flat truncation. In the test- 
body limit, v — 0, all the Hamiltonians coincide and for the equal-mass case, 
v = 1/4, their differences are largest. 



3-4 Functional representation of compact objects 

Before going on with the presentation of more dynamical expressions we 
will discuss in more detail the <5-function-source model we are employing. 
Although we are interested in both neutron stars and black holes, our matter 
model will be based on black holes because these are the simplest objects in 
general relativity and neutron stars resemble them very much as seen from 
outside. The simplest black holes are the isolated non-rotating ones. Their 
solution is the Schwarzschild metric which solves the Einstein field equations 
for all time. In isotropic coordinates, the Schwarzschild metric reads 

ds* = - (j^j ^ + (l + 4 , (69) 

where M is the gravitating mass of the black hole and r 2 = (a; 1 ) 2 + (a; 2 ) 2 + 
(x 3 ) 2 , dx 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 . It should be pointed out that the origin 
of the coordinate system r = is not located where the Schwarzschild singu- 
larity R = (radial Schwarzschild coordinate R) in Schwarzschild coordinates 
is located, rather it is located on the other side of the Einstein-Rosen bridge, 
at infinity. The relation between isotropic coordinates and Schwarzschild co- 
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ordinates reads R = r(l + ^§) 2 if MG/2c 2 < r and R' = r(l + §§r) 2 if 
< r < MG/2c 2 , where R' is another Schwarzschild radial coordinate ap- 
propriate for the geometry of the other side of the Einstein-Rosen bridge. 
The regimes < R < 2GM/c 2 and < R' < 2GM/c 2 are not accessible 
to isotropic coordinates. The harmonic radial coordinate, say here, g, relates 
to the Schwarzschild radial coordinate through R = g + MG/c 2 . Evidently, 
the origin of the harmonic coordinates is located at R — MG/c 2 which is a 
spacelike curve in the region between event horizon and Schwarzschild singu- 
larity. 

For two black holes, the metric for maximally sliced Brill-Lindquist initial 
time-symmetric data reads 

/ 1 _ PiG feG \ 2 / ^ ^, \ 4 

1 2^ \ 2 ,,2 , ( t , "1" , "2° \ , 2 



V + §& + ^) V 2r lC 2 2r 2 c 2 ' 



where the a a coefficients are given in Eq. (|4l)|) and where the f3 a coefficients 
can be found in [39] (notice d t r a = 0, initially). 

The total energy results from the ADM mass-energy expression 

Eadm = - jf ds ^ = J d 3 *AV = {ax + a 2 )c 2 , (71) 

where = 1 + „ aiG 2 + " 2G 2 and dsi = n l r 2 dQ is a two-dimensional surface- 
area element with unit radial vector n l = x l jr and solid angle element df2. 

Introducing the inversion map r[ = a 2 G 2 / Ac 4 ri or, rj = ria 2 G 2 /4c 4 r 2 
and ri = r^afG /4c r'i , where ri = x — xi, ri = |x — xi|, r[ = x' — xi, 
r[ = |x' — xi|, the three-metric at the throat of black hole 1, dl 2 = ^ 4 <ix 2 , 
transforms into 

dl 2 = W-d^ = ( 1 + 2% + 4 dx' 2 , (72) 

a 2 G 2 r' 

where r 2 — a \ c A + *12 with ri 2 = ri — r 2 . From the new metric function 

= 1 + 9 "/ G 2 + "^"rV 4 * ne P r0 P er mass of the throat 1 results in, taking 
into account v 2 = r' 1 a 2 G 2 /4c 4 r' 1 2 + ri2, 

(73) 



'20 

CtlO^G 

2n 2 c 2 



This construction as performed in Ref. [33j is a purely geometrical or vacuum 
one without touching singularities. Thus having the two individual masses 
mi and m 2 the gravitational interaction energy is obtained as E = Eadm — 
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(mi + Tn 2 )c 2 . Recall that this energy belongs to an initial value solution of 
the Einstein constraint equations with vanishing of both hj^ and particle 
and field momenta. In this initial conditions spurious gravitational waves are 
included. 

Let us introduce now point masses as sources for the Schwarzschild black 
hole. The stress-energy tensor density of test-mass point particles in a (d+1)- 
dimensional curved spacetime reads 

VOO = c 2 J2 m ^a,(tK(t)d a , (74) 

a 

where g^ v u a ^u av = g a ^ v u^u u a = -1, v£ = u£/u° a , and 5 a = <S(x - x a (t)) 
is the usual Dirac delta functions in d-dimensional flat space. In canonical 
framework, the momentum density Hi of the matter is given by 

7Ti = -if = cy^m a u ai 5 a . (75) 

a 

It is the source term in the momentum constraint. The important energy 
density in the Hamilton constraint reads 

- T*n» = ~c 2 m a u a ^S a = c 2 J2 m a u° a N8 a , (76) 

a a 

where n v is the timelike unit vector, n^n^ = —1 orthogonal to time slices t 
= const, n M = (— N, 0,0, 0); N is the lapse function, see [BJ. Independently 
from the form of the metric, the equations of motion are fulfilled and can be 
put into the form 

<V„u a „ = , (77) 

which means geodesic motion for each particle. 

The formal insertion of the stress-energy density into the Einstein field 
equations yields the following equations for the metric functions, for u ai = 0, 

_ . , 167rG \ -> 

- $A<f> = — g- ' ( 78 ) 



where 



lij =*^Sij, * = 1 + ^ d _\^ ■ ( 79 ) 



If the lapse function N is represented by 

X 



N = t , ( 8 °) 
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an equation for \ (be aware of the difference with \± m t ne Eqs. (|64[) and 
(j6"5")l ) results of the form, 

* A x = ^ rj — [ J2 m "X s °- ( 81 ) 

a 

With the aid of the relation, 

m 

it is easy to show that for 1 < d < 2 the equations for \P and \ do have 
well-defined solutions. Plugging in the ansatz, for d dimensions, 

G(d-2)P((d-2)/2) / ax a 2 A 
+ c 2 (d-l)7r( d - 2 )/ 2 ^r?~ a rg~V ' 1 J 

gives, for "mass-point 1" (mass-point seems the better notion compared with 
point particle or point mass because it is a fictitious particle only) , 

1 G(d-2)r((d-2)/2) / ai a 2 \\ . 

or, taking 1 < d < 2, and then taking the limit r\ — > 0, 

G(d-2)r((d-2)/2) « 2 \ A , . 



c 2 {d- l)7r( d - 2 )/ 2 



2 



Going over to (i = 3 by arguing that the solutions are analytic in d just results 
in the equation, cf. Eq. ([57]) , 

(86) 



" 1 + 

Tab 

where b ^ a and a, b = 1, 2, with the solution shown in Eq. (j4"6"|) . The ADM 
energy is again given by, in the limit d — 3, see Eq. (|71[) . 

Eadm = (ai + a 2 )c 2 ■ (87) 

Here we recognize the important property that although the Eqs. (|4"6"|) and 
(|4"T|) may describe close binary black holes with strongly deformed apparent 
horizons, the both black holes can still be generated by mass-points in con- 
formally related flat space. This is the justification for our particle model 
to be taken as model for orbiting black holes. We also will argue that bi- 
nary black holes generated by mass-points are orbiting black holes without 
spin, i.e. binary Schwarzschild-type black holes. We wish to point out that 
at the support of our ^-functions the physical spacetime is completely flat so 
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they can not be interpreted as local sources of gravity. They rather repre- 
sent wormhole geometries. The geometrical vacuum calculations in J33j are 
completely hnite, no infinite energies enter. The same holds with dimensional 
regularization where the formally infinite self-energies turn out to be zero. 
This is nicely consistent with the fact of the Dirac delta functions are living 
in flat space. 

Working in the sense of distributions, the dimensional regularization proce- 
dure preserves the important law of "tweedling of products" , [41] , F reg G re g = 
(FG) re g, and gives all integrals, particularly the inverse Laplacian, a unique 
definition. In the sense of analytic functions, all integrals are well-defined. A 
famous formula derived in |40j plays an all-over important role in PN calcu- 
lations, 



J d d x rfrf = 7T° 



rV a+d \ pfP+d\ r*f a+/3+d \ 
7d Y r <* r P = K_2_J±±_2_J±S 2 l r a+P+d 

r(-§)r(-f)r(^±f±^) 



It is well known that distributions or generalized functions like the Dirac 
delta function are boundary-value functions. To overcome distributional 
derivations like in 

aV -^ Pt ( ( ,_ 2) £i^)__^_ M , (89) 

where Pf denotes the Hadamard partie finie, it is very convenient to resort 
on the class of analytic functions introduced in [40] . 



r((d-e)/2) 



,e-d 



(90) 



e n d / 2 2*r(e/2) 
resulting in the Dirac delta function in the limit 

6 = lim e ^ <5 £ ■ (91) 
On this class of functions, the inverse Laplacian operates as 

-a-%= r[ }*- 2 r; )l2) y^=8^ (92) 

4^/22^(6/2 + 1) 

which is special case of the convolution property S e * S e > = <5 c + e ' which also 
results in the formula of Eq. (|88|) , and the second partial derivatives read 

d t d r^- d = Pf ( (d - 2 - e) (d-E)»y-frA . (93) 



No delta-function distributions are involved. Though the replacement in the 
stress-enery tensor density of 6 through S e does destroy the divergence free- 
ness of the stress-energy tensor and thus the integrability conditions of the 
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Einstein theory, the relaxed Einstein field equations (the ones which result 
after imposing coordinate conditions) do not force the stress-energy tensor to 
be divergence free and can thus be solved without problems. The solutions 
one gets do not fulfill the Einstein field equations but in the final limits of 
the e a going to zero the general coordinate covariance of the theory is recov- 
ered. This property, however, only holds if these limits were taken before the 
limit d — 3 is performed, |42j . For completeness we give here the terms which 
violate the contracted Bianchi identities, 



We wish to point out here a difference between ADM formalism and the 
harmonic coordinates approach. If in the harmonic coordinates approach the 
stress-energy tensor is not divergence free the relaxed field equations can 
be solved but the harmonic coordinate conditions will not be satisfied any 
further. This is different with the form we use the ADM formalism where the 
coordinate conditions are kept valid when solving the relaxed field equations. 
The relaxed field equations in the harmonic case include ten functions, which 
are just the metric coefficients, and the divergence freeness of the stress- 
energy tensor is achieved if on the solution space the harmonic coordinate 
conditions are imposed. In the ADM formalism in Routhian form the ten 
metric functions do fulfill the ADM coordinate conditions and equations of 
motion do follow from the Routhian. They, however, will not be the ones 
resulting from the Einstein field equations. Those will be obtained in the 
limits of e a — > only. 

The method of dimensional rcgularization has proven fully successful in 
both approaches, the Hamiltonian one and the one using the Einstein field 
equations in harmonic coordinates. However, another important difference be- 
tween both approaches should be mentioned. Whereas in the ADM approach 
all poles of the type l/(d — 3) cancel each other and no regularization con- 
stants are left, in the harmonic gauge approach poles survive with uncancelled 
constants, [S], [TO] - As found out, the difference is of gauge type only and can 
thus be eliminated by redefinition of the particle positions. On the other side, 
it shows that the positions of the mass-points in the Hamiltonian formalism 
are excellently chosen. Resorting to the maximally extended Schwarzschild 
metric, the spatial origin of the harmonic coordinates has Schwarzschild coor- 
dinate R = MG/c 2 inside horizon which can be reached by observers whereas 
the spatial origin of the ADM coordinates is located on a spacelike hyper- 
surface at R' = oo beyond horizon. The location of the origin of the ADM 
coordinates allows quite a nice control of the motion of the objects. 

The ADM coordinate system we are using in our article is called asymptot- 
ically maximal slicing because the trace of the extrinsic curvature of the t = 
const spacelike slices is not zero but decays as 1 /r 3 (in four spacetime dimen- 
sions) at spacelike infinity. It is closely related with the Dirac coordinate con- 




(94) 



a 
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ditions, (7 1 / 3 7 y ') j — 0, K\ = 0, which introduce maximal slicing. Recently, 
maximal slicing coordinates of the type introduced in Ref. [43] have proved 
useful in numerical relativity using moving punctures [44 . These coordinates 
are completely different from both the Dirac and ADM coordinates because 
those slices e.g. for Schwarschild black holes show crossings of the event hori- 
zon and settle down in the region between the Schwarzschild singularity and 
event horizon asymptotically. Only asymptotically these coordinates become 
rigidly connected to the Schwarzschild geometry. Nonetheless, moving punc- 
tures in numerical relativity are closely related with evolving Brill-Lindquist 
black holes. 



3. 5 PN expansion of the Routh functional 

In case of the full Einstein theory, a formal PN expansion of the Routh 
functional in powers of 1/c 2 is feasible. Using the definition hj^ = 16 ^ G hj^ , 
we may write 

oo ^ 

R [xl,p ai , hj?, dthj?) ~ c 2 J2m a = Y^ ^R 

a 71 — 

Furthermore, also the field equation for hjj can be put into a PN series form, 

(A - f ) h]? = £ ±D^[ x \ x t Pak , hl?,d t hl?) ■ 06) 
V c / n =0 

This equation has to be solved iteratively with the aid of retarded integrals 
which themselves have to be expanded in powers of 1/c. In higher orders, 
however, log-of-l/c terms will show up, [3J. 



.(95) 



3. 6 Near-zone energy loss versus far-zone energy flux 

The change in time of the matter Hamiltonian (it is minus the Lagrangian 
for the gravitational field) reads, assuming 1Z to be local in the gravitational 
field, 

dR or [dn- r on . ran- 
^=m=Jlm h + Jdvk Vh + JjK h > ^ 

where 

R = R[xi,p ai ,h,h] = JK(xi,p ai ,h,Vh,h) (98) 
with abbreviations 
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= / d 3 X, h = hj^ Vh = a fc/l TT j h = Q th TT (9g) 

Above, the equation for dR/dt is valid provided the equations of motion 

Pai = ~^ <= d^ (1 ° 0) 
hold. Furthermore, we have 

dK ■ fdK-- f f dK -\ d ffdK-\ , . 
;Vh+ / - r h= / V l—h)+- / l^-h) (101) 



dVh J dh J \dVh J dtj \dh 



dK 

\dVh )'" J dt \ dh 
Introducing the canonical field momentum 

c 3 _&R 

with abbreviation 7r = 7r^ T , and the Legendre transform 
„3 , c 3 



h. 



(102) 



the energy loss equation takes the form 



dH 
~dt 



f „ ( dK ,\ [ dK- f„ ( dK\ ■ f d I 'dK 



The application of the field equations 
dK 



h. (104) 



dh V {dVh)~dt{-BA)- ^ 
results in, employing the leading order quadratic field structure of K, 
dH f^fdKA r , dK ■ c 4 / ,„_,,,; 

^ = j v {wh h ) = v s Wh h= ^v s{vh)h (io6) 



dQr^hf , 



32irG Jf z 

where "fz" denotes the far zone (see e.g. section UTTj) . dfl the solid-angle 
element, and r the radial coordinate of the two-surface of integration with 
surface-area element ds = nr 2 df2. Here the further assumption has been 
made that the volume integrals in Eq. (|97[) may have the outer-most region 
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of the far zone as outer boundary. The expression 

'-sg/^W (107) 

is the well known total energy flux (luminosity C) of gravitational waves. 
The Newtonian and 1 PN wave generation fit into the above scheme of local 
Routhian density and far-zone as outer boundary which can be inferred from 



4 Binary point masses to higher PN order 

Most compact representations of dynamical systems are with Hamiltonians. 
Up to the 3.5 PN order, the Hamiltonian of binary point-mass systems is 
explicitly known, reading 

H(t) = rnxc 2 + m 2 c 2 + H N + -^H [1PN] + \h\ 2 pn] 

+ —H[ 2 .5PN](t) + —H[ 3PN ] + -fHl3.5pff](t) . (108) 

The non-autonomous dissipative Hamiltonians H^.spn] (t) and H^ 3 5PN j (t) 
are written as explicitly depending on time because they depend on the grav- 
itational field variables or, in case those are reduced to matter variables, on 
primed matter variables, see section l4~4l 

To simplify expressions like in section [5731 we go over to the center-of-mass 
frame Pi + P2 = and also define 

H = (H — mc 2 )//!, /i = vn\ui2 /w, m = mi + m 2 , v = fi/m, 

P = Pi/M> Pr = (n-p), q = (xi - x 2 )/Gm, n = q/|q| (109) 

with <v < 1/4 (y = test-body case, v = 1/4 equal-mass case) . 



4-1 Conservative Hamiltonians 

The conservative binary Hamiltonians read in reduced variables (the dissipa- 
tive Hamiltonians will be treated in section [4~4]) . see [47] . 

V 2 1 

H N = 7 L--, (110) 
2 q 
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H [1PN] = - l)p 4 - ^[(3 + ^)p 2 + + _L , (HI) 



-[(5 - 20i/ - 3t/ 2 )p 4 - 2^ 2 p 2 p 2 - 3^]- 
8 q 



+ 1[(5 + 8*V + 3^ 2 ]1 - 1(1 + 3z/)l , 1112) 



#[3P;v] - (-5 + 351/ - 70S + 35u 6 )p' 

+ -^[(-7 + 42i^ - 53j/ 2 - 5z/ 3 )p 6 + (2 - 3^)^ 2 p 2 p 4 

+ 3(l-^)^ 2 p 4 p 2 -5 l / 3 ^]i 

Q 

+ [^(-27 + 136i/ + 109i^ 2 )p 4 + ^(17 + 30^)^ 2 p 2 

+ 1( 5 + 43^)^ 4 ]1 
12 

25 ( 1 2 335 \ 23 2 \ 2 



-7T V V p 

8 \64 48 J 8 1 

85 3 2 7 \ 2l 1 .1 /109 21 2 \ . 1 . . 

tt 2 v vr>t\ — + - + 1" H-t(-113) 

16 64 4 / Fr q 3 L 8 V 12 32 / J ^ ; 

These Hamiltonians constitute an important element in the construction of 
templates for gravitational waves emitted from compact binaries. They serve 
also as basis of the effective one-body (EOB) approach, where with the aid 
of a canonical transformation the dynamics is put into test-body form of a 
deformed Schwarzschild metric, [51] . From the reduced Hamiltonians, where 
a factor of \jv is factorized out, the standard test-body dynamics is very 
easily obtained, simply by putting v = 0. 



4- 2 Dynamical invariants 

Dynamical invariants related to our previous dynamics are easily calculated 
within a Hamiltonian framework, |32j . Let us denote the radial action by 
i r (E,j) with E — H andp 2 = p 2 +j 2 /r 2 (p = p r e r +p v e v with orthonormal 
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i r (E,j) 



1 

2^ 



dr p r , 



(114) 



where the integration is originally defined from minimum to minimum radial 
distance. Thus all expressions derived hereof relate to orbits completed in this 
sense. From analytical mechanics it is known that the phase of the completed 
orbit revolution <P is given by 



<P d 

- = l + k = --i r (E, j ) 



and the orbital period P reads 

P 



d 



i r (E,j) . 



2-nGm dE 

Explicitly we get, for the periastron advance parameter k, 



1 3 



k — „.„<!+ „ 



c 2 j 2 



1 



1 



1 



7(7-2^)-^+ g (5- 21/) E 
1 , , E 



oiM^i + a -2( v )^ + a z{v) e z 



and for the orbital period, 
P 



2itGm (-2E) 3 / 2 
1 



1 f 1_1 

Efl 2 X 1 ~ 



{l5-v)E 



+ 



1 

c 6 



|( 5 - 2^) ( 2 ^ )3/2 - J-(35 + 30^ + 3^) E 2 
2 j 32 



(-2^ (-2g)Va 

a2<^) ^ oa 3 (v) : 

3 3 



■<u{v) E 3 



where 



ai(v) = 


5 /77 


2 V~2~ 


ai {v) = 


105 


- + 


asty) = 




a 4 (v) = 


lS< 21 



41 2 125 
64 ^ ~ ~3~ 

41 2 218 
64 71 " 3~ 

, i /1,,2\ 



v + -zr 
4 



45 



(115) 



(116) 



(117) 



(118) 

(119) 
(120) 
(121) 
(122) 
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These expressions have direct applications to binary pulsars, [45) . Explicit 
analytical orbit solutions of the conservative dynamics through 3 PN order 
are given in [46] . 



4-3 ISCO and the post- Newtonian framework 

The motion of a test-body in the Schwarzschild metric is known to have its 
innermost stable circular orbit (ISCO) at 6AlG/c 2 , in Schwarzschild coor- 
dinates. For test-bodies in rotating black holes (Kerr black holes) the ISCO 
lowers down up to MG/c 2 in case of direct motion and goes up to 9MG/c 2 
for retrograd motion, both motions in the equatorial plane. The ISCO of 
MG/c 2 is just the corotating case, where the Kerr black hole rotates as fast 
as the test-body is orbiting. In the both limiting cases of test-body motion 
MG/c 2 and 9MG/c 2 , the black holes have maximal spins. 

Within a Hamiltonian formalism the determination of the ISCO can be 
done straightforwardly. Just the following two equations have to be satisfied 
in the center-of-mass frame within the class of orbital circles (p r = 0), 

dH(r, J) = Q kal cirdes) d 2 H(r J) = Q 

or or A 

where r is the relative radial coordinate and J the orbital angular momentum. 
With the aid of the relation for the orbital frequency 

which holds for circular motion, the condition for the ISCO can also be put 
into the form 

— r-^=0. 125 

OLD 

Of course, one also could work with the Legendre transform F(u>) = H ( J) — 

Juj and put d — but the outcome would be the same (recall, dH = 
u>dJ, dF — —Jdio). Stability in the present approach with circular orbits 
means > or, ¥■ < or ^ < 0. 

In the context of approximation calculations, the main difference between 
the H(r, J) and H(u>) approaches is that r is not a coordinate invariant 
variable; in the case of approximately known H (r, J), this results in different 
values for the ISCO depending on the chosen coordinates. This is not the 
case with an approximately known H(u>) because it is coordinate invariant. 
Hereof, however, it does not follow that the ISCO calculated with H(uj) is 
more realistic than the ones calculated via H(r, J) rather the varieties of 
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ISCOs obtained via approximate H (r, J) show up the uncertainty in their 
true location. 

A test particle in the Schwarzschild spacetime on circular orbits has the 
reduced Hamiltonian 

1 3 2 27 , 675 4 3969 , ,,„„, 

= -2 a+ 8^ + 16^ + 128^ + l56^ + -"' (126) 

where 

v , , H(x)-mc 2 (GMuj\ 2/3 , . 

^ rnc? ' X= {—) • (m) 

The condition = yields i = l/6» 0.167 or, in Schwarzschild coor- 

dinates R = 6GM/c 2 . Evidently, ISCOs are located close to the reliability 
limit of PN expansions, [57] . 

Using the dynamical invariants, the angular frequency of circular motion 
can be written as 

1 + fc 

k-^circ — ^radial ^pcriastron — 27T — . (128) 

With the aid of the definition x = ( GM c " clrc ) , the binary dynamics yield, 
through 3PN order, i.e. c 2 H 3PN = H N + ^H [1PN] + ^H[ 2 pn] + ^H [3PN] , 

r'r / N X f3 I \ 2 (21 19 1 2 \ 3 

H 3PN (x) = h — I v \x + v H ir x 6 

6 yj 2 \8 24 J \16 16 48 J 

/675 / 34445 205 2 \ 155 , 35 ,\ 4/1 _, 
+ ( 128 + (-U52 + 192" J V+ + 10368* j ^ 

The ISCO, calculated with the aid of Eq. l(T2g|l . turns out to be x w 0.255, 
08], [37]. 

To likely improve the PN truncation a representation of the energy is 
helpful which in the test-particle limit is known to be a ratio of two simple 
polynomials e.g., 

1 — 4t 

e = (1 + H) 2 - 1 = - x Y~Tfa ( 13 °) 
For the binary 3 PN Hamiltonian, e turns out to be 

e(x) = -x[l- fl-^x- (3-^x 2 (131) 
n 5 /41 2 4309 \ 103 2 1 ,\ 

9 H — 7T 2 V + V v x A + 
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Applying to this expression the technic of padeing, i.e. putting it on a ratio- 
of-polynomials footing, the ISCO turned out to be, for equal-mass binaries, 
x w 0.198, [37]. 

To make contact to a recent discussion about the existence or non-existence 
of ISCOs for equal-mass binaries in a post-Newtonian setting [BJ we discuss 
the more general stability conditions for non-circular orbits using the Hamil- 
tonian H(r,p r , J), where J — p v . The crucial conditions for stability are 

Instability occurs if one of the >-signs turns to zero. This particularly means 
that the >-sign in (3] has to be zero. On the other side, in approximation 
calculations, where truncated series occur, the expression (3] can be zero 
without one of the expressions (1] or (2] being zero because the product of 
two series of the order n PN, which is of order 2n PN, is again truncated at 
n PN and thus can be zero without one or both of the factors being zero. In 
Ref. [5J the condition (3] has been given priority because it turned out to be 
coordinate invariant through 3 PN order (notice in this regard, p r — dW/dr, 
where W is an action). 



4-4 PN dissipative binary dynamics 

The leading order 2.5 PN dissipative binary orbital dynamics is described by 
the non-autonomous Hamiltonian, 49J, 

2G d 3 Qij(t) ( pupij p 2i P2] Gm 1 m 2 \ , , 

tt[2.5PN]{t) = ■=-g -T^ 1 , (133) 

5c° at 6 V mi m.2 ri2 / 



where 



(t) - m °( X '>a ~ \< 5 ^ ( 134 ) 



is the Newtonian mass-quadrupole tensor. Evidently, only after the Hamil- 
ton equations of motion are calculated the primed position and momentum 
variables resulting via Qij (t) from time differentiations and use of the equa- 
tions of motion are allowed to be identified with the unprimed position and 
momentum variables. The 3.5 PN Hamiltonian is known too, but it will not 
be given here because of quite lengthy expressions, [50]. Applications of the 
2.5 PN Hamiltonian can be found in, e.g. [S5], [S3], [5T], [55], where in Ref. 
[3"8] a transformation to the Burke-Thorne gauge (coordinate conditions) is 
performed. 
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Within the ADM formalism the action functional (i.e. the integral of the 
Lagrangian) of rotating bodies must have the following structure as long as 
the lengths of the spins are preserved in time, 

w = ( £ + £ s « ) n * ] + w- 1 d " x *ttW ( 135 ) 

- H ADM W a ,p ai ,S^\hJ^ 



TT 



Here, f}£ — fl a (i) — 5 e yfcAl(;)0')Ai(0(k)> Aj(i)(k)^a(j)(k) — 4x(k)(i) — 

Sij, tijk = (i - - k)(k - i)/2; and p ai , x l a , S$ = 5 a(i ), and 60^ = 
^CijkA a mrj\5A a mrty are the independent matter variables, where the index 

a again numerates the particles. Notice that O a are anholonomic variables 
related to the angle- type variables A a u^\ of the proper rotations (spins). 

Qa is the spin precession angular frequency vector of the ath particle. The 
equations of motion for the particles read, 

,i Sjdt'H ADM . ,,. Sjdt'H ADM , 100 , 

X ^ ) = 6 Pai (t) ' p «® = - 6*ffi (l36) 

m*) = SJ ?*™ M . ^°(*) = ^(*)^(t) , (137) 

5Sa (t) 

where the last equation results from the action functional through variation 

(i) 

with respect to O a ; for more details, see e.g. [M]. The Hamiltonian which 
generates both the evolution equations for all dynamical variables as well 
as the contraint equations through variation with respect to the Lagrange 
multipliers, the lapse and shift functions N and N l , respectively, is given by 

H = J d 3 x(NH- N l Hi) + E{yij], (138) 

where 



E hij] = j> dsihijj ~ Ijj.i) (139) 

is a surface integral at spacelike infinity i with dsi the two-dimensional 
surface-area element, [55], and 

n = n Rcld + n mattei , (uo) 

Hi = Hf cld + Hj natter (141) 
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with 

l^ wfield = _ 7 v 2R + J_ ^ _ ^ t (142) 



16ttG 



H 



field 



2c^ + tt*'^ 



(143) 



are the total Hamilton and (linear) momentum densities. After imposing the 
constraint equations H = Hi = and the coordinate conditions (|32|) - (|34|) . 
the energy expression S[7y] turns into the ADM-Hamiltonian, 



H 



r i n . c(i) l,TT ij 
• L aiPaii J a ) 'Hj ! /l TT 



(144) 



To linear order in the spin variables, the matter densities read (recall, 
5 a = <5(x — x a ); on simplicity reasons the index a will not show up in the 
following equations) 



-{np)cS- -t^T^.fe 



cpt 



mc — np 



(145) 



= PiS 



with 



and 



t 



l mk S lk 8 



PiPk 



np(mc — np) 



(j mk Sf + i mp 5 k l ) 1 ql S qp 8 



7 



np = -n^pp = (m 2 c 2 + "f lj PiPj) 
ki Si(iPj) 



1/2 



np 



-6 + 1 kl 1 mn 



S, 



m(iPj)PnPl ^ 

(np) 2 (mc — np) 



(146) 



(147) 



(148) 



where n M is the future directed unit vector field orthogonal to the t = const 
hypersurfaces, n^ = (— N, 0, 0, 0). The above expressions were shown to be 
correct up to (and including) the orders S/c 4 and S/c 2 in n mattcr and ?l™ attcr , 
respectively. 



Introducing a dreibein field e!j with e^'e 



Ijk and e (i ) fc e^ - 



a 



spin tensor SV fe 



(fe)(0 



'(fc) (0 ^ 



Sa can be introduced which fulfills the relation 



7 l V%Sfcz = 2%)%) = const, 



(149) 
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where 2Sf{\ = £ijkS(j)(k) ■ Crucial for our canonical formalism is the constancy 
in time of StqStiy Because of the symmetry property of the metric coefficients 
7y, the symmetric root of 7^ 

eaeij = 7^ , e»j = eji (150) 

can be taken for the dreibein held, i.e. &(j)k = &jk- 

To the order the formalism has been developed consistently, the following 
relations hold, 

Pat = f d 3 xK la " er , (151) 

JVa 

j aij = f d 3 x(.T i H; nattcr - x mr ttet ) = x\p ai - 4 Pai + s aim , (152) 

JVa 

where V a denotes the volume of particle a. Furthermore, 

i x a,Paj} = kj , {S a (i),Sa(j)} = <^i]kS a ( k ) , zero otherwise , (153) 
and the total linear and angular momenta respectively take the forms Pi — 
The crucial consistency relation reads, 

£ rrmattcr r 

-jyjr- = 2 N ^ T « (154) 

with 

^matter = j ^(^rnatter _ ^matter) _ ^ggj 

It is fulhlled to the needed order, 

VjT i:j = -^-5 + t^ k + 0(G). (156) 

In asymptotically flat spacetimes the Poincare group is a global symme- 
try group. Its generators P p and J pv are conserved and fulfill the Poincare 
algebra, see e.g. [55] . 

{P",P u } = 0, (157) 
{P p , J pa } = -Tf p P a + rf a P p , (158) 
{J pv ,P pa } = —if p J pa + ri pp r a + rf p J pv - if p J pp . (159) 

The meaning of the components are energy P° — H/c, linear momentum 
P l = Pi, angular momentum J y = Jy, and Lorentz boost J' l0 /c = K l — 
G % — tP l . A center-of-mass vector can be defined by X 1 = c 2 G l /H. This 
vector, however is not a canonical position vector, see e.g. [54] . The energy 
H and the center-of-mass vector G l = Gi have the representations 
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Gi = -T^G S^^ = -^j/ d an^ X %- 5i ^ (161) 

where i° denotes spacelike infinity, r 2 dQw is the two-dimensional surface- 
area element, and n the radial unit vector. The two quantities, H, G l , are 
the most involved ones of those entering the Poincare algebra. 

In terms of three-dimensional quantities, the Poincare algebra reads, see 
e.g. [35], with Jij = e i3k Jk, 

{P i ,H} = {J i ,H} = 0, (162) 

{■h, Pj} — £ijk Pk 7 {Ji, Jj} — £ijk Jk > (163) 

{J i ,G j } = s ijk G k , (164) 
{Gi,H} = Pi, (165) 

{G i) P j } = ±H8 iJ , (166) 

{Gi,G 3 } - -^e ijk J k . (167) 

The Poincare algebra has been extensively used in the calculations of PN 
Hamiltonians for spinning binaries, |18j . |19j . Hereby the most important 
equation was (|165p which tells that the total linear momentum has to be a 
total time derivative. Once this equation has even fixed the kinetic ambiguity 
in non-dimensional regularization calculations, [56] . The kinetic ambiguity 
got also fixed by a Lorentzian version of the Hadamard regularization based 
on the Fock-de Donder approach, [57] , 



5.1 Approximate Hamiltonians for spinning binaries 

All the Hamiltonians, and the center-of-mass vectors too, given in this section 
have been derived or rederived in recent papers by the author and his collab- 
orators employing general relativity in canonical form, [17] - |21j . adjusting 
it to the motion of binary black holes. 

The Hamiltonian of leading-order (LO) spin-orbit coupling reads 




(168) 



and the one of leading-order spin(l)-spin(2) coupling is given by 
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ttLO 

U S 1 S 2 



EE 

a b^a 



2c 2 rl 



[3(S a .n ab )(S b .n Qb )-(S Q -S fc )] 



(169) 



where r ab n ab = x a — x b , a ^ b and a,b — 1,2. The more complicated Hamil- 
tonian is the one with spin-squared terms because it relates to the rotational 
deformation of spinning black holes. To leading order, say for spin(l), it reads 
(details of derivation are given later) 



H s° = ^iSr [3(Si • n ia )(Si • n 12 ) - (S x • Si)] 



(170) 



The LO spin-orbit and spin(a)-spin(6) center-of-mass vectors take the form 



(171) 



Within the conservative 3 PN dynamics for spinless point masses, the center- 
of-mass vector has been calculated in [SS] . Applications of the LO spin Hamil- 
tonians can be found in e.g. [58], [59], [60]. Other references treating the LO 
spin dynamics are e.g. [61], [62], [63], [64]. For applications of the next-to- 
leading-order (NLO) spin dynamics, presented straight below, see, e.g, [16] . 
[65]. 

The Hamiltonian of the NLO spin-orbit coupling reads, r = ria, 



SO 



= -G 



((pi x Si) • m 2 ) 



5m 2 pf 3(pi ■ p 2 ) 

8mf 4m| 



3pj 

4mim 2 



3(pi • n i2 )(p 2 • m 2 ) 3(p 2 • ni 2 ) 2 



G 



G 



((p 2 x Si) • ni 2 ) 



((Pi x Si) • p 2 ) 



2mirn 2 

(Pi-P2) . 3(pi • ni 2 )(p 2 • m 2 ) 



TOim 2 mim 2 
2(p 2 • m 2 ) 3(pi • ni 2 



G 



2 ((pi x Si) -n 12 ) 



mim 2 



llm 2 5m 2 



G 



2 ((P2 x Si) ■ n 12 ) 



2 

6mi 



mi 
15to 2 



(1 



2) 



and the one of NLO spin(l)-spin(2) coupling is given by 
G 



H 



NLO 
S1S2 



2mim 2 c 4 r 3 



[6((p 2 x Si) • ni 2 )((pi x S 2 ) • ni 2 ) 



(172) 
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3 

+ -((pi x Si) • n 12 )((p 2 x S 2 ) • n i2 ) 

- 15(Si • ni 2 )(S 2 • ni 2 )(pi ■ ni 2 )(p 2 • n 12 ) 

- 3(Si • n i2 )(S 2 • ni 2 )(pi ■ p 2 ) + 3(Si • p 2 )(S 2 • ni 2 )(pi • n i2 ) 

+ 3(S 2 • pi)(Si • nia)(p2 • n i2 ) + 3(Si ■ Pi)(S 2 • ni 2 )(p 2 • m 2 ) 

+ 3(S 2 • p 2 )(Si • ni 2 )(pi • n i2 ) - 3(Si ■ S 2 )(pi ■ ni 2 )(p 2 • n i2 ) 

+ (Si • Pl )(S a • p 2 ) - i(Si • p 2 )(S 2 ■ Pl ) + i(Si ■ S 2 )( Pl ■ p 2 )] 
3 

+ 9 o [~((pi x Si) • ni 2 )((pi x S 2 ) • n 12 ) 
Lm\T a 

+ (Si • S 2 )( P i • ni 2 ) 2 - (Si • n 12 )(S a • Pi)(pi • n i2 )] 
3 

+ „ 2 J ~((P2 x s 2) • n i2 )((p 2 x Si) • n 12 ) 

+ (Si • S 2 )(p 2 • m 2 ) 2 - (S a • ni 2 )(Si • p 2 )(p 2 • n 12 )] 

+ 6(mi t7 )G2 [(Si • S a ) - 2(Si • n 12 )(Sa • n 12 )]. (173) 

The calculation of the LO and NLO order S l2 -Hamiltonians needs more in- 
formation about the source terms than given in Eqs. (|145p and (|146j) . To 
achieve the LO + NLO Sf-Hamiltonians, the following additional source in 
the Hamilton constraint is needed, 

^static = + ^rnn^ l ^ nk JuAklSl{m 

where ; i and , i denote three-dimensional covariant and partial derivatives, 
respectively, and where 

Q ij = Y k l 3l l mn Si km S lnl + ^S 2 7 y ' , (175) 
2S 2 = j lk ^ l SnjS lkl = const. 

Qi is the quadupole tensor of the black hole with number 1 resulting from its 
rotational deformation. Herewith, beyond the previously shown LO Hamil- 
tonian (|170p . the NLO Hamiltonian comes out in the form, employing the 
Poincare algebra for unique fixation of all coefficients, 



Post-Newtonian methods: Analytic results on the binary problem 



33 



H 



NLO 



G 



1712 ( C \2 



4mim 2 



3m 2 



3m 2 



^l(pi-n) (Sx-n) (p^Si) 
3 



4mim2 



n 2 S 2 



/ \2 2 2 / Q \* 

9 



4mim2 



p^(Si-n)" 



4m 2 

3 
rn 



(pi • p a ) S 2 - —j ( Pl • p 2 ) (Si • n) 2 - — , ( Pl • n) (p 2 • Si) (Si • n) 



4m 



2m\ 



2 (p 2 • n) ( Pl • Si) (Si ■ n) + — ( Pl • n) (p 2 • n) S 2 



'i 
15 



4m 2 



G 2 m 2 
2c 4 r 4 



pi ■ n) (p 2 • n) (Si • n) 

14m 2 



9(Si • n) 2 - 5S 2 + ^^(Si • n) 2 - ^S 2 



771 1 



mi 



(176) 



The spin precession equations of the Hamiltonians HgQ° and H^ 2 LO have 
been calculated also in the papers [23] and [21], respectively, where the first 
paper [23] has benefited from paper [20]. The final spin precession equation 
of the second paper [24], Eq. ([60]) , deviates from the corresponding one in 
[2T] . A detailed inspection has shown that the last term in Eq. (|62[) of [25] has 
wrong sign, [26j . Using the correct sign, after redefinition of the spin variable, 
agreement with the Hamiltonian of Ref. [2T] is achieved. 

The NLO order spin-orbit and spin(a)-spin(6) center-of-mass vectors take 
the form 



G 



NLO 
SO 



E^L(p° xS «) 



/-i NLO 



iTVLO 



+ EE 

a b=jta 

+ EE 

a b=£a 



mi,G 
Ac 4 m a r ab 



((Pa x S Q ) • n afe ) 5Xa +Xb - 5(p Q x S Q ) 



Tab 



G 

C^ab 



3 1 

-{p b x S a ) - ^{n ab x S Q )(p 6 • n a6 ) 



((p a X S a ) • tl ab ) 



(177) 



2^4 



2m 2 G 
c 4 mi 



3 (Si -n 12 ) 



(xi + x 2 ) + —i- (3xi - 5x 2 ) - 

12 ° r 12 



(178) 
(Si -ni 2 )Si 

(179) 
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summing up through 2 PN order in the spin parts, we obtain, also see |56j . 

Grt i f~< , /-i l r i piO i riNLO i r<NLO i r<NLO i pJVLO 

— Lrjv + ljripjv + '^2PA r + '-J3PA r + ^'S'0 + lj 'SO + tj SiS , 2 +tj S 2 +"5? 

1 . (2l80) 

Numerically, spins of black holes can be counted of order 1/c (maximum), 
thus the spinless parts are taken up to the 3 PN order. 

The currently known conservative binary Hamiltonians for spinning black 
holes through order 1 /c can be summarized as follows 

H = H n + Hipn + H 2 pn + H3PN 

+ w so + H s 1 s 2 + H sf + H sl 

, ttNLO , ttNLO , ttNLO , tjNLO 
+ ^SO + W SiS 2 + ^6\ 2 + 

+ #piSf + -^W? + ^piS? + #p a Sj 

+ HpiStSl + H p 2 S 2 Sj + H Pl S 2 Sf + H P2SxSi 

+ H S 2 S 2 + H SlS 3 + H S2S 3 . (181) 

The Hamiltonians -ffg4 and iJ S 4 in the approximation in question turned out 
to be zero. 



6 Lorentz-covariant approach and PN expansions 

The Lorentz-covariant approach has found a quite thorough presentation in 
[2J. So we will not go into so many details as in the canonical approach 
presented in the previous sections. 

The Einstein field equations are given by 

G A {g K \id a g K \,d a d[3g K x) = ^ , (182) 

where g K \ and are the four-metric and the stress-energy tensor of the 
matter (e.g., fluid), respectively. The contracted Bianchi identities yield the 
four-dimensional equations of motion (EOM) for the matter, 

V V G» V = -► V„T^ = (EOM), (183) 

where V„ denotes the four-dimensional covariant derivative. The Landau- 
Lifshitz form of the Einstein field equations fits very well into Lorentz- 
covariant schemes. It takes the form, [3J, 

d x d K U^(g a p) = ^-rH{g a0 , d 7 g afj ) , (184) 
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with U^ vXk = g^g AK - g^ A g !/K and g^" = y/^g^", where g denotes the 
determinant of the metric tensor, t^l is known as the Landau-Lifshitz stress- 
energy (or energy-momentum) pseudo-tensor of the gravitational field. It is 
unique in the sense of symmetry and dependence on the metric coefficients 
and its first derivatives only. The equations of motion now read 

d u d x d K U^ = -» d v T»l = (EOM) (185) 

with 



r, 



= -gT^ + —^t^(g aP ,d 7 g^) . (186) 

Applying the condition of harmonic coordinates 

d v ^ v = or, d„H^ = with H* v = sf^gg^ - rf v , (187) 

where rf v denotes the Minkowski metric, the Einstein field equations in 
Landau-Lifshitz form read (so-called relaxed field equations because they 
do not imply the equations of motion; rather the condition of harmonic co- 
ordinates implies the equations of motion), 

1 £\ f ' 

rf^d a dpH' tu = — j- t" u , (188) 



where 



and 



-gT^ + (189) 



= -H^dadpH^ + d a H^dpH ya + ^ v g all d x H aT d T H px (190) 

+ \{2g^g^ - g^g^)(2g XT g pa - g Xa g Tp )d a H x ° d H T ? 

- g^gprdxH^BaH^ - g ua gp T d x W d a H^ + g afi g XT d x m a d T H^ . 

The A^ v object starts with quadratic non-linearities of the gravitational field. 
It is another stress-energy pseudo-tensor of the gravitational field. 



6.1 PM and PN expansions 

The formal retarded solution (resulting under the condition of no incoming 
radiation) of the inhomogeneous wave equation (1188)) reads 

H^(xi,t) = -^- f rfVr^V, t - |X - ; c 2 )|x - x'l' 1 . (191) 
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A post-Minkowskian (PM) expansion in powers of G can now be introduced 
in the form. 



H^(x,t) = Y / G n H^ ] (x,t;c). (192) 



n=l 

If additionally the virial theorem holds, 

GM V 2 



Rc 2 



,:2 



(193) 



where respectively M, R, and V are typical masses, radii, and velocities of 
the system in question, the PM expansion may be further expanded into a 
PN series in powers of 1/c. Yet, because of the retardation structure of the 
solution, a PN expansion is achieveable only in the near and far zones and 
this even only in a generalized form with log-c terms showing up at higher 
orders starting from 4 PN, i.e. (1/c 2 ) 4 = 1/c 8 , on. 

Let us assume now that the matter source is bounded by a sphere with 
radius R centered in the origin of the coordinate system and that for the 
typical gravitational wave length A the relation A >> R holds. The near 
zone is then defined by |x| << A. The formal PN expansion (near-zone PN 
expansion) is defined by 

AC _ 00 _ ( —'\\ n r f) n 
^(x,*) = -^5:^/rf 3 x'-^^(x' ) t;c 2 )|x-xr- 1 . (194) 

n— 

Additionally, the expansion 



1 

i 

n=-l 



rHx',; ;c 2 )= E ^sTw^*) ( 195 ) 



applies. In the far zone, where r = |x| >> A holds, a formal PN expansion 
(far-zone PN expansion) yields 

C(M) = ~ / d 3 x'-^r-(x',t-^; C 2 )(x'.„r. (196) 

n— 

Here, the expansion 

oo 1 



n=-l C 



applies. The Eqs. (filMf and (fT95|) on the one side and (fT96| and (fl"97f on 

the other are somewhat simplified in the sense that they do not show up 
log-c terms at higher orders in the expansions which result from the badly 
defined integrals of non-compact support. On the other side, as the expres- 
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sions stand, they are mathematically not defined at all. Details can be found 
in the works by Blanchet, Damour, Will, and collaborators; particularly see 
the contribution by L. Blanchet in this volume; for tail terms, also see our 
section [ 



6.2 PN expansion in the near zone 



Up to the 2 PN order the metric coefficients read 

2 „ 2 



gm = _ 1+ V _ V <+ {x + ViVi + 



9oi — — oVi r-Rz 



4 - 



(198) 
(199) 
(200) 



With the following choice of the matter variables, respectively mass, mass- 
current, and stress density, 



^00 _j_ rpii 



Cf C 

the 2 PN potentials can be put into the form 

rf 3 



= T ij 



F(x, i) = G rct {-4ttGct} = G [ j^-M^ t - |x - z|/c) 

J l x - z l 



(201) 



(202) 



Vi = G xet {~4irGa z } , 



(203) 



Wij = G rct {-47rG((Ty - 5ij(x k k) - diVdjV} 



(204) 



Ri = G rct <j -^G{Va l - V^) - 2d k VdiV k - -d t Vd t V } , 



(205) 



X = G rct -inGVau + 2V i d t d i V + Vd 2 t V 



(206) 
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The potentials of the orders 2.5 PN and 3.5 PN are radiation-reaction 
potentials. They are most compactly given under Burke-Thorne coordinate 
conditions, reading 



C/ reac (x,f) 



189 



G 

„5 



21 



- x 



(207) 
(208) 

(209) 

(210) 
(211) 



The used definitions read yV = y l y\ y<«> ee tfi = STF(y^), and e.g., M]l 
indicates the seventh time derivative of M^. Explicitly, the 1 PN metric 
including the gravitational radiation reaction through 3.5 PN order is given 

by 



where the source multipole moments are given by 



Sij = / d 3 y e km<l y 3>k (Jr, 



goo = -1 + + [/ reac ) + 1 [9 t 2 X - 2(7 2 - 4[/[/ r ' 

sM = -^i+ur c ), 



where the potentials have the integral representations 

d 3 y 



U(x,t) =G J 
Ui{-K,t) = G 



|x-y| 

d 3 y 
|x-y| 



*(y,t)> 



0i(y,t), 



x(x,t) = G /" d 3 y |x-y|a(y,i). 



(212) 
(213) 
(214) 

(215) 
(216) 
(217) 



These integrals are well defined. Evidently, multipole expansion and PN ex- 
pansion nicely fit together; see also [56]. 
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6.3 PN expansion in the far zone 
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In the far zone, the multipole expansion of the transverse-traceless (TT) part 
of the gravitational field, obtained by algebraic projection with Pjjfc m (n), 
reads, e.g. [67] , 



^fzTT( X ' ^) 



G P i]km (n) 



Eli 



1=2 



i M m (t--)N- 

^ m kmi 3 ...ii\ l c I iV *i 



+ (?) " (7^^ S lL3...,( i -7)^^3...,| ! (218) 



where the leading mass-quadrupole tensor takes the form, e.g. 



Mi 



(-7) 



2GV 



/•OO 




11" 


/ 


In ( ^ 




/o 


V26/ 


+ 12 



M^'(*--f-«) + - (219) 



with 



2Gm, 
r* = r H tt- m 



showing a leading-order tail term. Notice the modification of the standard PN 
expansion through tail terms. The Eq. (|218p nicely shows that also multipole 
expansions in the far zone do induce PN expansions. 

The gravitational luminosity is generally given by {H^ TT = —hJ Tiz ), 



(220) 



Through 1.5 PN order, the luminosity explicitly reads, 
G 



C(t) 



y - 

5c 5 ^ \ c 2 

n=0 v 



C n {t) 



G_ 

5? 



M [3l M [31 



1 



l^c[3]o[3] 



(221) 



On reasons of energy balance, for any representation of the Einstein theory, 
the time-averaged energy loss has to fulfill a relation of the form 



- < 



d£(t 



dt 



> = < C(t) > 



(222) 
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where the time averaging procedure takes into account typical periods of the 
system. The derivation of this equation in section 13.61 is known to be valid 
for the first two radiation emission and reaction levels. 



7 Energy loss and gravitational wave emission 

The energy flux to n PN order in the far zone, denoted n PN(fz), implies 
energy loss to (n+5/2) PN order in the near zone, denoted (n+5/2) PN(nz). 
Hereof it follows that energy-loss calculations are quite efficient via energy- 
flux calculations. Because of this we will apply the balance property between 
emitted and lost energies to some PN orders to easily derive the energy loss 
from the energy flux. In general, only after averaging over orbital periods the 
both expressions will coincide (see Eq. (|222p V In the case of circular orbits, 
however, this averaging procedure is not needed. 



7.1 Orbital decay to 4 PN order 

The binding energy of our binary system on circular orbits is given by (iE c i r 
Therefore, for the energy loss to 4 PN order, we get 



dE r 



/'- 



dt 



C 



32c E 
~5G 



v 2 x 5 



1247 35 



V 336 12 



v x + Anx 



3/2 



(223) 



where the 1.5 PN(fz) energy flux is taken from Ref. [69] where also the 2 
PN(fz) energy flux can be found; for the 3.5 PN(fz) energy flux see [70], [7Tj . 

Taking into account the Eq. (| 1 29[1 we obtain a differential equation for x 
which is easily solved with accuracy 1/c 8 . In terms of the dimensionless time 
variable 



> = rr=-(tc-i), 



5Gm 



(224) 



where t c denotes the coalescence time, the solution reads [69] 



X = -T- 1 / 4 



1+f 



743 



V 4032 48 



(225) 



Taking into account the relation between phase and frequency ^ = -ir- 
respectively 7p = — -x 3 / 2 , the phase evolution results in 



I T 5/8 



^ + \ -1/4 _ 3^ T -3/8' 

8064 96 7 4 



(226) 
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7.2 Gravitational waveform to 1.5 PN order 



The radiation field can be decomposed into two orthogonal polarization 
states. The polarization states h+ and h x are defined by 



h+ = ^(uiUj - ViVj)hJ^ ', (227) 

h x = -(nvj + v iUj )hjy, (228) 

where u and v denote two vectors in the polarization plane forming an orthog- 
onal right-handed triad with the direction n from the source to the detector. 
The detector is directly sensitive to a linear combination of the polarization 
waveforms h + and h x , namely 

h(t) = F+h+(t)+F x h x (t), (229) 

where F + and F x are the so-called beam-pattern functions of the detector 
depending on two angles giving the direction — n of the source as seen from 
the detector and a polarization angle specifying the orientation of the vectors 
u and v around that direction. 

For our binary system, the two polarizations h+ and h x are chosen such 
that the polarization vectors u and v lie respectively along the major and 
minor axis of the projection onto the plane of the sky of the circular orbit, 
with u oriented toward the ascending node, the point at which black hole 1 
crosses the plane of the sky moving toward the observer. The result, to 1.5 
PN(fz) order, reads [75j (the 2 PN(fz) wave form is given therein too) 



2G(j,x 
' l +,x — — 5 



(230) 



where the plus polarization is given by 

Hf = -(l + c?)cos2^, (231) 
H [i/2] = _£i^[(5 + c 2) cos ^ _ 9(1 + c£)cos3# (232) 



8 m 

H [ l ] = g[19+ 19c- - 2c\ - i/(19 - 11c? - 6cf)]cos2^> 

- |sf(l + cf)(l-3^)cos4^, (233) 

4 3/21 = 4^ — {[57 + 60c? -4- 2,(49 - 12c? - c \ )]cos^ 
192 m 

27 

- — [73 + 40c? - 9cf - 2^(25 - 8c? - 9cf)]cos3^ 
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625 

+ — (1 - 2i/)s?(l + c?)cos5^} - 2tt(1 + c 2 )cos2^ , (234) 
and the cross polarization by 

H [ ° ] = -2c 4 sin2V> , (235) 
H l * /2] = —Sid — [sin^ - 3sin3V>] , (236) 

-^x 1 = §[ 17 - 4c ? - ^( 13 - 12c?)]sin2^ - !|c 4 s 2 (l - 3j/)sin4^, (237) 

4 3/2I = ^^{[63-5cf-2K23-c?)] S i^ 
96 m 

27 

- y [67 - 15c, 2 - 2i/(19 - 15c 2 )]sin3V> 
625 

+ — (1 - 2z/)s 2 sin5^} - 47r Ci sin2-0 , (238) 

where Cj = cosi and = sim and i denotes the inclination angle between 
the direction of the detector, as seen from the binary's center-of-mass, and 
the normal to the orbital plane which is assumed to be right-handed with 
respect to the sense of motion so that < i < ir. dm = mi — m2, and the 
phase variable -0 is given by 

ip = 4> - 3x 3/2 ln (j^J , (239) 

where 4* is the actual orbital phase of the binary, namely the angle oriented 
in the sense of motion between the ascending node and the direction of black 
hole 1 (<f> = mod 2tt when the two black holes lie along u, with black hole 
1 at the ascending node). The logarithmic phase modulation originates from 
the propagation of tails in the wave zone. The constant scale xo can be chosen 
arbitrarily; it relates to the arbitrary constant b in the Eq. (|219p . For details 
on higher order PN levels, see e.g. [2]. 
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